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O , Abstract. We show that the ground state of a Bose-Einstein condensate of atoms 

with hyperfine spin f — 2 can be either spin aligned, condensed into pairs of atoms 
coupled to = 0, or condensed into triplets of atoms coupled to F = 0. The complete 
^ I phase diagram is constructed for f — 2 and the generic properties of the phase diagram 

are obtained for / > 2. 



Q ; PACS numbers: 03.75.Fi, 03.65.Fd 

If atoms in a Bose-Einstein condensate (BEC) are trapped by optical means [1], 
their hyperfine spins (or spins) are not frozen in one particular direction but are 
essentially free but for their mutual interactions. As a result, the atoms do not behave 
as scalar particles but each of the components of the spin is involved in the formation 
■ of the BEC. This raises interesting questions concerning the structure of the condensate 
and how it depends on the spin exchange interactions between the atoms. 

Such questions were addressed in a series of theoretical papers by Ho and co- 
workers [2] who obtained solutions based on a generating function method. In the 
case of spin-1 atoms the problem of quantum spin mixing was analyzed by Law et al. [3] 
who proposed an elegant solution based on algebraic methods. It is the purpose of 
^ ■ this paper to point out that a wide class of many-body hamiltonians appropriate for 
^ the problem of interacting bosons with spin can be solved through algebraic techniques 

which have found fruitful applications in nuclear physics [1] as well as in other fields 
of physics (see, e.g. Ref. [5]). The main result derived in this paper is that an exact 
solution is available for spin values / = 1 and / = 2 (for any value of the number of 
atoms A^) which allows the analytic determination of the structure of the ground state 
of the condensate. For spin values / > 2 solvable classes of hamiltonians give insights 
into the generic properties of the phase diagram. 

We consider a one-component dilute gas of trapped bosonic atoms with arbitrary 
(integer) hyperfine spin /. In second quantization the hamiltonian of this system has 
a one-body and a two-body piece that can be written as (we follow the notation of 
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Ref. [3]) 

where ^ = 1, Ma is the mass of the atom, and and are the atomic field 
annihilation and creation operators associated with atoms in the hyperfine state |/m) 
with m = — /, . . . , +/, the possible values of all summation indices in ([1]). The trapping 
potential Vtrap is assumed to be the same for all 2/ + 1 components. The coefficients 
^rnxm^mzm^ follow from the interaction between atoms which is assumed to be of short- 
range, two-body character, 

?7(f„ X,) = 6{x, - X,) E ^fI/'; FM){f- FM\, (2) 

FM 

where |/^;FM) is the combined state of the atoms i and j with total spin F, and 
Up = A-nfi^ap / with ap being the s-wave scattering length in the F channel. The 
assumption underpinning the form ([2]) is rotational invariance of the hamiltonian in 
hyperfine-spin space. 

We assume in this paper that the scattering lengths in the different F channels 
are comparable and that, in first approximation, the interaction strength between the 
bosons is independent of F. In that case the dominant part of the hamiltonian ([1]) is of 
the form 

n, = ni + x' I ^l.^U^m.^m.d'x, (3) 

mim2 

and is symmetric under any interchange of the spin-component indices. Under this 
assumption the condensate wave functions for each spin component (f)m{x) {m = 
— /, . . . , -|-/) can be approximated by a single wave function which satisfies the 
Gross-Pitaevskii equation associated with the dominant hamiltonian [3]. Furthermore, 
the atomic field creation and annihilation operators at zero temperature can be 
approximated by 

^l-bl(f){x), ^^^bm(j){x), m = -f,...,+f, (4) 

where 6^ and b}^ are annihilation and creation operators associated with the entire 
condensate, satisfying the usual boson commutation rules 

In this approximation the entire hamiltonian ([1]) can be rewritten as 

H = eb^ ■b + -J2Mb^ X b^]^^^ ■ [by<b]^^\ (6) 
2 p 

where the coefficients e and up are related to those in the original hamiltonian through 
integration over x, viz. up = u'p J \(j){x)\'^d^x. The notation x in Eq. implies the 
coupling to a given spin F and projection M, 

[b^ X 6t]CJ) = Y: {fmfm'\FM)blbl„ (7) 
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where (■ ■ ■ ■ | ■ ■) is a Clebsch-Gordan coefficient Furthermore, the dot ■ denotes a 
scalar product, 

.f^^ (-)^V2F + l[f ^ X t^]S,°\ (8) 

for tensor operators of rank F. The definition of the adjoint operator bm = 
{—)^^^b_^m ensures that b^ is an annihilation operator with transformation properties 
under rotations that are the same as those for the creation operator b\^ [7J. With the 
above definitions we have that b'^ - b = X^m Hn^m is the number operator which counts 
the total number of atoms in the condensate. 

To derive the solvability properties of the hamiltonian ([6]), we ffist determine its 
algebraic structure by introducing the bilinear operators bj^bm'- From Eq. ([5]) one finds 
the commutation relations 

which can be identified as those of the unitary (Lie) algebra U(2/ + 1) [7]. Exactly 
solvable hamiltonians with rotational or SO (3) invariance are now found by the 
determination of all Lie algebras G satisfying U(2/ + 1) D G D S0(3). The canonical 
reduction of U(2/ + 1) is of the form 

U(2/ + l) D S0(2/ + l) D S0(3). (10) 

[For / = 3 there is an additional exceptional G2 algebra between S0(2/ + 1) and 
S0(3) which for the symmetric representations of U(2/ + 1) considered here does not 
add anything to the discussion.] The relevance of a chain of nested algebras of the 
type (|TOl) is that it defines a set of commuting operators and with it a class of solvable 
hamiltonians. Consider in particular the hamiltonian 

H' = aiC'i[U(2/ + 1)] + a2C2[V{2f + 1)] 

+ 6C'2[SO(2/ + l)] + cC'2[SO(3)], (11) 

where ai, 02, b, and c are numerical coefficients and CnfG] is the n^'^-order Casimir 
operator of the algebra G which satisfies the property that it commutes with all 
generators of G [8]. Solvability of the hamiltonian (fTTl) follows from the fact that it 
is written as a sum of commuting operators, a property which indeed is valid for the 
Casimir operators associated to any chain of nested algebras such as (ITO!) . The Casimir 
operators appearing in Eq. (ITTll are known in closed form, 

(7i[U(2/ + l)] =iV, 
C'2[U(2/ + 1)] =N{N + 2f), 

C'2[SO(2/ + 1)] = - (2/ + l)f ° ■ r° + N{N + 2/ - 1), 

C2[SO(3)] = [iFiF + 1) - /(/ + 1)] f^-f^ + fif + (12) 

F 

in terms of the operators a/ = [^^ ^ ^^]m'' ^^^1 T_f m = ^ ^]m^- Equations (IT^ show 
that the solvable hamiltonian f|TT]) is a special case of the general hamiltonian ([6]) with 



Bose-Einstein condensates of atoms with arbitrary spin 



4 



coefficients e and up that are Hnear combinations of ai, 02, b, and c according to 
e =ai + (2/ + l)a2 + 2/6 + /(/ + l)c, 
i/i. = 2a2 + 26 + + 1) - 2/(/ + l)]c, F 0, 

uo =2a2-4/6-2/(/ + l)c, (13) 
Tlie eigenvalues of the hamiltonian (fTTi) are 

E'{N,v,F) = aiN + a2N{N + 2f) + bv{v + 2f - 1) + cF(F + 1). (14) 

The allowed values of v are f = A^, — 2, . . . , 1 or 0, as can be obtained from the 
U(2/ + 1) D S0(2/ + 1) branching rule [8]. The quantum number v corresponds 
to the number of bosons not in pairs of bosons coupled to F = 0, and is known 
as seniority P, [10]. The allowed values of the total spin F are obtained from 
the S0(2/ + 1) D S0(3) branching rule which is rather complicated but known in 
general [11]. The / = 2 example is discussed below. 

The generic solvability properties of the original hamiltonian ([6]) now follow from a 
simple counting argument. For atoms with spin / = 1 the solvable hamiltonian ([TTl) has 
three coefficients ai, 02, and c [since SO(2/+l)=SO(3)] while the general hamiltonian ([6]) 
also contains three coefficients e, z/q, and 1^2. [Note that the coupling of two spins to 
odd F is not allowed in the approximation (|4]) of a common spatial wave function, so no 
z/i term occurs.] For atoms with spin / = 2 both the solvable and general hamiltonian 
contains four coefficients (ai, 02, b, and c versus e, z/q, z/2, and z/4) which can be put into 
one-to-one correspondence. Hence the general hamiltonian ([6]) is solvable for / = 2. 
The same counting argument shows that it is no longer solvable for / > 2. 

The case of interacting / = 1 atoms was discussed by Law et al. [3J who identified 
the existence of two possible condensate ground states: one with all atoms aligned to 
maximum spin F = N and a second with pairs of atoms coupled to F = 0. Whether 
the condensate is aligned or paired depends on a single interaction parameter which in 
our notation is c. With the technique explained above we can also derive the phase 
diagram for atoms with spin / = 2. The results are exact and valid for arbitrary 
A^. The entire spectrum is determined by the eigenvalue expression ( JT^ together with 
the necessary branching rules. In particular, the allowed values of total spin F for 
a given seniority v are derived from the SO (5) D SO (3) branching rule [1] given by 
F = 2r, 2r - 2, 2r - 3, . . . , r + 1, r with t = v,v - 3,v - 6, . . . and r > 0. 

It is now possible to determine all possible ground-state configurations of the 
condensate. This problem has been considered in the study of the spectral features 
of quantal systems with random interactions [12]. We note that the character of the 
ground state does not depend on the coefficients since the first two terms in the 
expression ([T^ give a constant contribution to the energy of all states. Although this 
contribution is dominant according to our earlier assumptions, the spectrum generating 
perturbation of the hamiltonian is confined to the last two terms and depends solely on 
the coefficients b and c which are related to the original interactions up according to 

6 = ;^(-7z/o + lOz/2 - 3Z/4), c = -^(-z/2 + z/4), (15) 
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Figure 1. Diagrams of tlie different phases of a Bose-Einstein condensate of atoms 
with spin / = 2 characterized by a ground state (vq, Fq) where vq is the seniority of the 
ground state and Fq is its total spin. The total number of atoms N is even in (a) and 
odd in (b). The grey area corresponds to a ground state with {vq, Fq) = {N — 3 + 6,0) 
which only occurs for (5 = ±1 and disappears in the limit iV — > cx). 



The following exact finite- results are found where the ground state of the condensate 
is characterized by a seniority Vq and a total spin Fq. 

(i) is even. We introduce N = 6k + 26 with k integer and 6 = —1,0,-1-1. The 
possible ground-state configurations have {vo,Fq) = (0,0), {N,2N), {N,2\6\), or 
(A^ — 3 -|- 5, 0), the latter existing only for 6 = ±1, 

(ii) is odd. We introduce N = Qk + 3 + 26 with k integer and 6 = — 1, 0, -|-1. The 
possible ground-state configurations have (fo, Fq) = (1, 2), (3, 0), (A^, 2A^), (A^, 2\6\), 
or (A^ — 3 -|- (5, 0), the latter existing only for 6 = ±1. 

The phase diagram displays a richer structure than in the / = 1 case as is shown in 
Fig. [H We observe first of all the presence of the aligned phase where the seniority is 
maximal, Vq = N, and all spins are aligned, Fq = 2N. Secondly, we have a low-seniority 
(paired) and consequently low-spin phase. For even A^ this corresponds necessarily to 
{vo,Fo) = (0,0). For odd A^ there must be at least one unpaired atom leading to the 
ground-state configuration {vo,Fq) = (1,2); alternatively, however, it might consist of 
a triplet of atoms which is coupled to total spin Fq = leading to the ground-state 
configuration {vo,Fq) = (3,0). The (1,2) and (3,0) phases are divided by the line 
b = 3c/ 7. The paired and aligned phases are separated by the line 



2N{2N + 1) 



{2N -2){2N + 3) 



(16) 



N{N + 3) ' (A^-l)(A^ + 4) ' 

for A^ even or odd respectively, which in both cases tends to b = —4c for A^ ^ oo. 

So far we have recovered the aligned and paired phases also encountered for 
interacting / = 1 atoms (although the paired phase is somewhat more intricate for 
/ = 2 due to the possible presence of a triplet of atoms coupled to F = 0). For / = 2 
a third phase occurs for negative b and positive c characterized by high seniority {i.e. 
unpaired) and low total spin, {vo,Fo) = {N,2\6\). Finally, for 6 = ±1 there exists a 
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pathological region in the phase diagram characterized by {vq, Fq) = {N — 3 + S,0) (see 
Fig. [1]). It is separated from the high-seniority, low-spin region by the line 

Ji(* + 3)| 

M2N + S)' * ' 

which tends to 6 = for — oo. Hence this region disappears in the large- limit. 

We conclude that the ground state of a BEC consisting of atoms with spin / = 2 
can be of three different types: (i) a maximum- seniority spin-aligned, (ii) a low-seniority 
low-spin, or (iii) a maximum- seniority low-spin configuration. Note that 'seniority' in 
this context refers to number of atoms that are not in pairs coupled to F = 0. 

Since the hamiltonian (ITT]) is solvable for / = 2, all eigenstates, and in particular 
the three different ground states, can determined analytically. The general expressions 
given by Chacon et al. [13| reduce to 

N/2 



\v 


= N,F 


= M 


= 2N) 


oc 


\v 


= 0,F-- 


= M -- 


= 0) 


oc 


\v 


= N,F 


= M 


= 0) 


oc 



dt.dt)-/^|0) 



atxat](2).at)'^^'|0), (18) 

where the / = 2 atoms are denoted as d bosons. In the second of these expressions 
it is assumed that is even and in the third that = 3k; other cases are obtained 
by adding a single boson or an F = pair. The are the so-called traceless boson 
operators [13] which are defined as (see also Chapt. 8 of Ref. [H]) 

«L = dl - ^^yqr^^m- (19) 



We emphasize that (1181) are the exact finite- A^ expressions for the eigenstates of the 
hamiltonian (ITTl) . Since in the large- A^ limit the traceless boson operators a}^ become 
identical to dl^, we arrive at a simple interpretation of the three types of configurations: 
(i) spin-aligned, (ii) condensed into pairs of atoms coupled to F = 0, and (iii) condensed 
into triplets of atoms coupled to F = 0. 

How will these features evolve with increasing spin / of the atoms? For arbitrary 
interaction strengths up in the different F channels the hamiltonian is not solvable. 
By imposing / — 2 conditions on up it can be brought into the form (fTTl) and this gives 
an idea of the structure of the general phase diagram by constructing a two-dimensional 
slice of it. For example, for atoms with spin / = 3 the elimination of ai, 02, b, and c 
from Eq. ( fT3|) yields the condition llz/2 — 181/4 + Tuq = 0. For / > 3 more conditions on 
Up are found. If all conditions are satisfied, the phase diagram in b and c with 

-7z/o + lOz/2 - 3z/4 1 , . 

b= : — ^ , C= — (—1^2 + ^4:), 20 

14(2/ + 1) ' 14^ ^ ^ ^ 

has properties similar to those in the / = 2 case. The analysis requires the knowledge 
of the multiplicity d[J\F)^ {i.e., the number of spin-/ atom states with seniority v 
coupled to total spin F) which can be derived from the SO (2/ + 1) D SO (3) branching 
rule [11]. We find that for sufficiently large even A^ there are four competing ground 
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states with {vo,Fq) = {N,fN), {N,0), (0,0), and (2,2), the latter of which disappears 
as a ground state in the large- iV limit. For sufficiently large odd the four competing 
ground states have {vq, Fq) = {N, fN), {N, 0), (1, /), and (3, f2 = f mod 2), the latter 
two being separated by the line b = [/(/ + 1) — /2(/2 + l)]c/(4/ + 6). The results 
correspond to what is found in the / = 2 case and lead to an essentially identical (6, c) 
phase diagram. 

Finally, we point out that the appearance of exact seniority ground states requires 
weaker conditions on up than those that have been discussed so far. In fact, the spin- 
aligned configuration (A^, fN) is always an eigenstate of the general hamiltonian ([6]) 
because the F = fN state is unique. Furthermore, it can be shown [TO] that seniority 
is a good quantum number if the interaction strengths up satisfy [//3J conditions 
only (where \_x\ is the largest integer smaller than or equal to x). For all cases 
of any conceivable interest for BECs, this reduces to no condition on the strengths 
Up for / = 1,2 or just a single one for / = 3,4,5. So there is at most a single 
condition required for all eigenstates to carry exact seniority and for the results of 
this paper to be valid. Nevertheless, the determination of the complete phase diagram 
for / > 2 with unconstrained interaction strengths up remains a problem worthy of 
further investigation. 
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